Magnetoasymmetric current fluctuations of single-electron tunneling 
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We determine the shot noise asymmetry of a quantum dot under reversal of an external magnetic 
field. The dot is coupled to edge states which invert their chirality when the field is reversed, leading 
to a magnetoasymmetric electrochemical potential in the nanostructure. Surprisingly, we find an 
exact relation between the magnetoasymmetries corresponding to the nonfinear conductance and 
the shot noise to leading order in the applied bias, implying a higher-order fluctuation-dissipation 
relationship. Our calculations also show a magnetoasymmetry of the full probability distribution of 
the transferred charge. 



PACS numbers: 73.23.-b, 73.50.Fq, 73.63.Kv 



I. INTRODUCTION 

Recently, it has been theoretically 

demonstratedii^i^ii^i^i^ and experimentally 

verified^i^'^'iSiiii'i^ii^ that the Onsager-Casimir reci- 
procity relationsi^ii^ cannot be, in general, extended to 
mesoscopic transport far from equilibrium. Departures 
of microreversibility at equilibrium have been related to 
the interaction of the nanostructure with an environ- 
ment driven out of equilibrium.^^ In both cases, it is 
found that the current flowing through the system is 
not invariant under reversal of the external magnetic 
field, leading to the observation of magnetoasymmetries 
solely due to the asymmetric properties of the internal 
electrostatic potential. Thus, the effect is induced purely 
by electron interactions. 

Previous literature has discussed the size of the mag- 
netoasymmetries for the electric conductance. Here, 
we are concerned with the shot noise, which is known 
to offer a complementary and, quite often, unique 
tool to probe electronic transport in quantum corre- 
lated nanostructuresiii Shot noise has been shown to 
reveal electronic entanglement detection in mesoscopic 
interferometerSfiSi dynamical spin blockade in dots at- 
tached to ferromagnetic leadsfi^ quantum shuttling in na- 
noelectromechanical systemSf^S, nonequilibrium lifetime 
broadenings in cotunneling currents, ^"'^ and quantum co- 
herent coupling in double quantum dots,^^ just to men- 
tion a few. 

It is worth noting that the linear conductance and the 
equilibrium noise are related each other via the linear 
fluctuation-dissipation theorem, which is a general state- 
ment about the response of a system near equilibrium and 
its dynamical fluctuations induced by random forces. ^'^ 
For mesoscopic conductors the theorem is expressed as 
the Johnson-Nyquist formula between equilibrium cur- 
rent fluctuations and the linear conductance.^^ Both the 
linear conductance and the current fluctuations at equi- 
librium obey reciprocity relationsi^i At nonequilibrium, 
however, Onsager's microreversibility is generally not ful- 
filled. Notably, our calculations explicitly show an exact 
relation between both magnetoasymmetries correspond- 
ing to the noise susceptibility and the nonlinear conduc- 



tance to leading order in the applied voltage bias, im- 
plying a higher-order fluctuation-dissipation relationship. 
Recent work a^^'^^ find similar fluctuation relations. 



II. THE SYSTEM 

We show our result considering a simple but paradig- 
matic mesoscopic system — a quantum dot in the 
Coulomb blockade regime for which the charge is quan- 
tized and transport is blocked at low temperature unless 
charging energy is supplied by external voltage. We 
study a dot coupled to two chiral states^i^ (filling factor 

V — 1) propagating along the opposite edges of a quan- 
tum Hall conductor, as shown in Fig. [TJ For positive 
magnetic fields B > carriers in the upper (lower) edge 
state move from the left (right) terminal to the right (left) 
terminal. The current flow is reversed for B < 0. Cou- 
pling between the dot and the edge states takes place via 
tunnel couplings S71 and ^172 and capacitive couplings 
Ci and C2. 

We consider a single energy level Eq which can be ex- 
ternally tuned with a gate voltage. {Eq is a kinetic en- 
ergy invariant under B reversal). The dot lies deep in 
the Coulomb-blockade regime for which Eq + e^/C (with 
C = Ci -I- C2) is well above the electrochemical poten- 
tials /ic( = Ep + eVa {a = L, R) of both left (L) and 
right (i?) leads, Ep being the common Fermi energy. We 
assume that each edge state is in equilibrium with its 
corresponding injecting reservoir. Therefore, they act ef- 
fectively as massive electrodes with well defined electro- 
chemical potentials. Without loss of generality, we take 

VI = —Vr — V/2 with V the applied bias voltage. 

For temperatures ksT and voltages low enough, 
ksT^eV e^/C (but fc^T sufficiently high to neglect 
Kondo correlations), the charging energy e^/C is the 
dominant energy scale of the problem and double occu- 
pation in the dot is negligible. In the master equation 
approach, the charging state of the dot is given by an 
integer number of electrons, n. We assume that quan- 
tum coherence between states with differing n is lost, as 
occurs in the pure Coulomb blockade regime. For defi- 
niteness, we consider spin polarized carriers, although the 



2 





FIG. 1: Sketch of a quantum dot coupled to chiral edge states 
via tunnel and capacitive couplings. The sample is driven out 
of equilibrium with an external bias Vl = —Vr = V/2. 



model can be easily extended to the spinful case. Thus, 
the dot can be either in the empty state (n = 0) or in the 
occupied state (n = 1) with instantaneous occupation 
probabihty Pn{t) at time t. The main transport mech- 
anism is via single-electron tunneling (71,72 ^ ksT/K) 
and we can write quantum rate equations for po{t) and 



Pi 



-r+ r- 
r+ -r- 



Po 

pi 



(1) 



In a more compact form, one has p — Mp, where p is 
the vector with components p = {po Pif" and the 2x2 
matrix M can be obtained from Eq. ^ . We note that the 
columns of M add to zero to fulfill the condition po +Pi = 
1 at every instant of t. The elements of M represent 
transition rates, F', to tunnel on (e = -t-1) and off the 
dot (e = — 1). The total rates are given by F^ — F^. 
For S > we find (see Fig. [U F+ = 7i/(+S,Vl), 
r+ = I2f{+B, Vr), F- = 7i[l - /(+B, Vl)] and F^ = 
72(1 — /(+^, Vr)] where the occupation factors are given 
by /(B, 14) = 1/(1 + exp (^d(S) - eVadlkBT) for Ep - 
0. The electrochemical potential in the dot, Hd{B), is self- 
consistently found from the electrostatic configuration, 
which depends on the B direction,-'' For B > 0, one has 
^Jid{+B) = Eo- eCiVL/2C - eC^VRl^C. 

For S < the chirality of the edge states is in- 
verted and the nonequilibrium (Vl 7^ Vr) polariza- 
tion charge changes accordingly. Injected electrons in 
the upper (lower) edge state are now predominantly 
emitted from lead R (L). Therefore, the rates now 
read FJ 



function^ 

Sc^piu;) = 2 J e^"*[(/„(i)/^(0)) - dt , (2) 

where (/„) is the averaged current across the a junction. 
Its stationary value can be found from the steady state 
solution of Eq. ([T]), which reads po — F^/F and pi = 
F+/F. As a consequence. 



{la) = e[F+po - r„pi] = 



F+F" - F"F+ 



(3) 



where F = F+ -^F^. 

The correlator {Ia{t)Ip{{))) in Eq. ^ can be obtained 
from the conditional probability (t) that the state n 
is occupied at time t > Q when the dot was in the state 
n' aX t = Q. Within our scheme the quantum regression 
theorem holds, and Pn,n' (t) obeys the same equations as 
Pn{t). Hence, the eigenvalues of M completely determine 
the dynamical behavior of p. This can be seen from the 
noise expression Sap = S^^Sap + 5'^^, with the 
Schottky noise produced by correlated tunneling through 
a single junction and 



Slp{^) 



JGw fa;), (4) 



where the "Green function" matrix for Eq. ([T]) reads, 

G(w) = -Re(ia; + M)-i. (5) 

Care must be taken with the w = limit since M is 
singular. The matrix M has two eigenvalues, namely, 
Ai = with an eigenvector given by the stationary so- 
lution p = (F-/F F+/F)^ and A2 = -F, which de- 
scribes a charge excitation in the system. Therefore, 
we can now use in Eq. ([5]) the spectral decomposition 
M — X]i=i 2 -^jUEjU"^, where is a matrix with the 
element (i, i) equal to 1 and the other elements are ze- 
roes, the i-th column of U being the i-th eigenvector of 
M. The Ai contribution cancels out the term {Ia){Ip)- 
As a result, G{lu) = -A2UE2U-7(^^^ + ^2) and we find 
the shot noise S = —Slr. at = 0, 



2? 



n^R 



F3 



(6) 



l2f{-B,VL), FJ. = ^-^f[-B,VR), F^ = For B > and large bias such that eV > Eq^UbT, we 



72 [1 - f{-B,VL)] and F^ = - f{-B,VR)], where 
^id{-B) ^Eo- eC2VL/2C - eCiVR/2C. We stress that 
the rates depend, quite generally, on the B direction and 
differ for Ci 7^ C2. As demonstrated in Ref. [s^, the 
magnetoasymmetry of the polarization charge leads to a 
magnetoasymmetric addition energy of the quantum dot 
and, as a result, the current average is an uneven func- 
tion of B. We can thus ask ourselves whether the current 
fluctuations also exhibit this effect. 

Current fluctuations are characterized by the power 
spectrum, Sap{uj), of the current-current correlation 



have F+ « 71, F~ w 0, F+ « and F^ w 72. Substi- 
tuting these values in Eq. ([5]) we recover the double bar- 
rier case for noninteracting electrons, S = 26^7172(7^ -I- 
tD/t^: where 7 = 71 + 72. 



III. RESULTS 

Figure [5] shows results for the averaged current / — 
(II) = —{Ir) and the Fano factor F = S/2eI as a func- 
tion of V. The current is exponentially suppressed at 
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FIG. 2: (Color online) Current (black curves) and Fano fac- 
tor (red curves) as a function of bias voltage for 71 = 72 = 
0.005Eo/h, kgT = 0.1_Bo, Ci = 0.6 and C2 = 0.4. 



low V, increases at ~ CEq/bCi and reaches the hmit 
value 67172/7 for large voltages. When B is reversed, 
the current now increases at ~ CEq/ eC2- Since we 
choose Ci > C2 the /(V, —B) is shifted to larger voltages 
compared to I{V,+B). As a consequence, the differen- 
tial conductance is generally i?-asymmetric. The Fano 
factor is Poissonian at small V since transport is dom- 
inated by thermal activated tunneling. For increasing 
V , the noise becomes sub-Poissonian and reaches satu- 
ration for large V . The crossover step from Poissonian 
noise to sub-Poissonian noise has a width which depends 
on ksT. Like the current, the crossover center shifts to 
larger voltages when B is reversed, thus yielding a mag- 
netoasymmetric Fano factor. 

Since the Fano factor depends on both the noise and 
the current and these are asymmetric under B reversal, 
we plot in Fig. [3] the magnetic field asymmetry for the 
noise alone, which we define as <i>s — [S{+B) — S{—B)\/2. 
The asymmetry vanishes for small V , fulfilling the On- 
sager symmetry. At large V the noise saturation value 
is independent of the B direction since this limit cor- 
responds to noninteracting fermions. As a result, the 
asymmetry vanishes. The asymmetry becomes maxi- 
mal for intermediate voltages. Importantly, the asym- 
metry increases with the capacitance asymmetry, rj — 
(Ci — C2)/C. Therefore, the current fluctuations are 
magnetoasymmetric only in the case where the electro- 
static coupling of the dot with the edge states leads to 
an asymmetric screening of charges. 

Furthermore, we consider the case of low voltages, 
eV <C fc^T. Then, we can expand / and S in powers 
of V, 



/ = GiV + G2V^ + OiV"^) , 
S = So + SiV + 0{V^). 



(7) 
(8) 



We have checked that both the linear conductance Gi 
and the equilibrium noise Sq are even functions of B. 
They satisfy the linear fluctuation dissipation theorem, 
'S'o = ^ksTGi. While Gi depends on the equilib- 
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FIG. 3: (Color online) Shot noise magnetoasymmetry for var- 
ious capacitance asymmetries 77 = (Ci — C2)/(Ci + C2). We 
take 71 = 72 = 0.005So/fi and fcsT = Q.lEo- Inset: Mag- 
netoasymmetries of the leading-order nonlinear conductance 
and nonequilibrium noise at V = O.OSiSo /e. 



rium potential in the dot, the nonlinear conductance 
term G2 is, in general, a function of the screening 
electrostatic potential4=- This potential need not be 
B-symmetric— Thus, the magnetoasymmetry = 
[G2{B) — G2(— J5)]/2 acquires a finite value: 



6^7172 ?7_ 

8 7 (fcsT)2 



rsech 



^0 
2kBT 



tanh 



2kBT 



(9) 



For high fc^r, Eq. ^ yields — since thermal 

fluctuations are B-symmetric. determines the rate 

at which the current magnetoasymmetry increases with 
voltage, thus showing that magnetoasymmetries are a 
truly nonequilibrium effect. 

The nonequilibrium noise at linear response [i.e., 5*1 in 
Eq. ([5])] is also B-asymmetric. To leading order in V, we 
find 



$c, = AkBT<^( 



(10) 



where $5^ = [Si{B) - Si{-B)]/2. This is a relevant re- 
sult of our work. Remarkably, we obtain the same func- 
tional dependence for the magnetoasymmetries of noise 
and current in the leading-order nonlinearities. We nu- 
merically confirm this prediction for a small value of V 
(see inset of Fig. [S]). This result can be related to the 
nonlinear fluctuation-dissipation theorem, which reads 
Si = AkBTG2- It has been derived in Ref. '36' for meso- 
scopic conductors at arbitrary voltages within the frame- 
work of full counting statistics. Thus, it is shown''^ that 
the nonlinear fluctuation-dissipation theorem holds even 
for interacting particles assuming time-reversibility (no 
magnetic fields). The nontrivial difference, however, is 
that in our theory microreversibility is broken due to the 
combined effect of magnetic fields and interactions'^ but 
still Eq. PH]) holds. While detailed balance conditions 
are shown to hold far from equilibrium in the absence of 
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FIG. 4: (Color online) Charge probability distribution mag- 
netoasymmetry for 71 = 72 = 0.05Eo/h, ksT = 0.4iJo, 
eV = — 4.5-Bo, and various values of ry. Inset: P{N) for two 
field orientations and rj = 0.4. 

magnetic fields,— i^i^^ the same relations can not, gener- 
ally, be established when microreversibility is brokeni^^ 
Despite this, we obtain an unexpected symmetry rela- 
tion between the conductance and noise response mag- 
netoasymmetries. 

We note in passing that Eq. ^TU\i is not a general- 
ized of the fluctuation-dissipation theorem for which non- 
linear fluctuations and the response are related via a 
(nontrivial) effective temperature, as in glassy systems;^ 
since the prefactor for both linear and nonlinear the- 
orems is the same. Our result also differs from more 
general fluctuation theorems obeyed by full probability 
distributions .— 

To analyze higher-order terms one should consider fluc- 
tuations of the screening potential, which are beyond 
the scope of a mean-field approximation. Nevertheless, 
for classical Coulomb blockade effects the local poten- 
tial fully screens the excess charges and quantum fluctua- 
tions are absent. Therefore, our model system is perfectly 
suitable for further extensions. We note that the Fano 
factor magnetoasymmetry, = [F{B) — F(—B)]/2, is 
quadratic in voltage: 

^^^H^J^I^Vsech^-^tanh-^. (11) 
72 '\2kBTj 2kBT 2kBT ^ ' 

A complete characterization of current fluctuations is 



given by the full counting statistics,^ which yields the 
entire probability distribution P{N) of the transferred 
charge during the measurement time t^. We follow the 
method of Bagrets and Nazarov^ to assess the cumulant 
generating function S{x) = —^n^j^P{N)e''^^. With- 
out loss of generality, we count charges in the R lead. 
Therefore, we make the substitutions — + Fj^e*-^ and 
— > r^e~*-^ in the off-diagonal elements of M and 
S{x) is derived from S{x) = io'^2(x)- We calculate 
P{N) within the saddle-point approximation, valid in the 
limit Iq — > C!C^2 and determine the magnetoasymmetry of 
P{N). We show $p = [P(iV, +B) - P{N, -B)]/2 in Fig. 
m |$p| increases for increasing capacitance asymmetry 
and vanishes around N ~ 7to/4. This point corresponds 
to the mean current / = 67/4 for a dot symmetrically 
coupled (71 =72) in the limit eV ^ ksT. 



IV. CONCLUSIONS 

To summarize, we have investigated magnetoasymmet- 
ric current fluctuations of a Coulomb-blockaded quantum 
dot. It is well established that the linear fluctuation- 
dissipation theorem makes an equivalence between linear- 
response functions to small perturbations and correla- 
tion functions describing fluctuations due to electric mo- 
tion. In this work, we have found a similar fluctuation- 
dissipation relation that predicts an exact equivalence 
between the leading-order rectification and noise magne- 
toasymmetries, valid in the presence of external magnetic 
fields. Such relation has been very recently shown to de- 
rive from fundamental principles. Moreover, we have 
shown that the full probability distribution associated to 
the How of charges is, generally, magnetic- field asymmet- 
ric. Since the effect studied here relies purely on inter- 
action, it should be observable in many other systems 
exhibiting strong charging effects. 
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